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Infrared behaviour and fixed points in Landau gauge QCD 
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We investigate the infrared behaviour of gluon and ghost propagators in Landau gauge QCD by 
means of an exact renormalisation group equation. We explain how, in general, the infrared momen- 
tum structure of Green functions can be extracted within this approach. An optimisation procedure 
is devised to remove residual regulator dependences. In Landau gauge QCD this framework is used 
to determine the infrared leading terms of the propagators. The results support the Kugo-Ojima 
confinement scenario. Possible extensions are discussed. 
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Charged particles have long-range nature, whereas lo- 
calised objects are neutral. For local quantum field 
theories one derives more generally that every gauge- 
invariant localised state is singlet with respect to the 
unbroken global charges of the gauge symmetry. Thus, 
in QCD without Higgs mechanism any localised physi- 
cal state must be colourless. This extends to all phys- 
ical states only if a mass gap is present. Then, colour- 
electric charge superselection sectors cannot arise, and 
every gauge-invariant state must also be a colour singlet. 
These are the signatures of confinement. 

Within covariant linear gauges, the necessary condi- 
tions for confinement were formulated more than twenty 
years ago by Kugo and Ojima They express the ne- 
cessity of a mass gap and the conditions for avoiding 
the Higgs mechanism in a covariant continuum formula- 
tion of QCD in terms of local fields, i.e., for the local 
measurement of (colour) quantum numbers. Both con- 
ditions are encoded in the critical infrared exponents of 
gluonic and ghost correlations in Landau gauge QCD 0- 
For the mass gap, the massless transverse gluon states of 
perturbation theory must be screened nonperturbatively, 
leading to an infrared suppression for the gluon propaga- 
tor. In turn, the global gauge charges entail the infrared 
enhancement of the ghost correlations leading to their 
infrared dominance. 

The infrared behaviour just described has been found 
in solutions of truncated Dyson-Schwinger equations Q . 
Since then, the same qualitative behaviour has emerged 
in a variety of nonperturbative approaches including an 
increasing number of lattice simulations |J| and investi- 
gations based on Fokker-Planck type diffusion equation 
of stochastic quantisation [5(. Quantitatively, the lat- 
ter leads to the same infrared exponents as the Dyson- 
Schwinger equations, which reflects an equivalence be- 
tween the two formulations under certain conditions 
Interesting in their own right, such interrelations 
are especially useful for finding ways to go beyond the 
particular approximations employed in each case. 

Here we present an investigation of the infrared sector 
of Landau gauge QCD by means of a nonperturbative 
flow equation 0, Q . Central to the approach is the ef- 



fective action T^, where quantum fluctuations with mo- 
menta p 2 > k 2 are already integrated out, and k denotes 
an infrared cut-off scale. The variation with k leads to a 
flow equation for T^, thereby interpolating between the 
classical action in the ultraviolet and the full quantum 
effective action in the infrared where the cut-off is re- 
moved. Moreover, it has been shown that truncations 
can be systematically improved by means of an optimi- 
sation procedure The formalism has its particular 
merits in gauge theories as it is amiable to approxima- 
tions being non-local in momenta and fields, precisely the 
behaviour we expect in QCD. So far it has been applied 
to Landau gauge QCD for a determination of the heavy 
quark effective potential [Tol ] and effective quark interac- 
tions above the confinement scale ; for a brief review 
and further applications in Yang-Mills theories see 01 
and respectively. 

Conceptually, the flow comprises a system of coupled 
differential equations for all Green functions which only 
involves dressed vertices and propagators. Hence, the 
correct RG-scaling is displayed by the full flow and trun- 
cations with the correct symmetry properties. This is 
one of the main contradistinctions to Dyson-Schwinger 
equations or stochastic quantisation which involve both 
bare and dressed Green functions. This makes it often 
difficult in the latter two approaches to implement the 
correct RG-scaling behaviour within truncations. 

In this letter we explain how the full infrared momen- 
tum behaviour of Green functions is extracted within a 
fixed point regime. In Landau gauge QCD, this allows for 
a simple determination of infrared coefficients for gluon 
and ghost propagator as well as the running coupling. 
Results for general regulators are presented and inter- 
preted in the context of an optimisation of the flow. The 
systematic extension of the truncation is discussed. 

We begin by introducing the classical gauge fixed ac- 
tion for QCD in a general covariant gauge in four dimen- 
sions. Including also the ghost action, it is given by 



S cl = i / tr^ 2 + i / {D^f + I C- d^D^ ■ C. (1) 



The flow equation approach relies on a momentum cut- 



off for the propagating degrees of freedom. Here, the 
momentum regularisation for both, the gluon and the 
ghost, is achieved by S c i — > iS c i + ASk with 

AS k = | J A£ Rf v Al + Jc a R ab C b . (2) 

The functions R implement the infrared momentum 
cut-off at the momentum scale p 2 w k 2 . The scale- 
dependence of the regularisation in (J2J) induces a flow 
equation for the effective action. With cf> — (A, C, C) 
and t — Ink, the flow equation takes the explicit form 

dtT k [cb] = \Jd A P G%[4>](p,p) d t R%{p) 

- Jd 4 P G ab m P ,p)d t R ba (p), (3) 

where G[cj)}(p,q) = (T k 2) + Ry\p,q) denotes the full 
regularised propagator. The flow is finite in both, the 
infrared and the ultraviolet, by construction. Effectively, 
the momentum integration in @ only receives contribu- 
tions for momenta in the vicinity of p 2 m k 2 . Conse- 
quently it has a remarkable numerical stability. The flow 
solely depends on dressed vertices and propagators lead- 
ing to consistent RG-scaling on either side of ©. The 
flow equations for propagators and vertices are obtained 
from @ by functional derivatives. 

To verify the Kugo-Ojima confinement criterion we 
have to determine the momentum behaviour of the prop- 
agators in the deep infrared region 

k 2 « p 2 « A QCD . (4) 

In this domain, quantum fluctuations are already inte- 
grated out and the physical quantities as well as gen- 
eral vertex functions will show no k dependence. More- 
over, @ entails that we are already deep in the infrared 
regime. No new physics arises and the flow equation © 
only constitutes the shifting of the transition regime at 
momenta about k 2 towards smaller scales. Accordingly, 
Green functions 1^ = S^T^/S^i ■ ■ ■ 5<fi n can be param- 
eterised as 

lf>(p? « Aq GD ) = z n t^( Pl ,p 2 /k 2 ) , (5) 

where pi stands for the momenta p\, ...,p n - Here z n are 
possibly fc-dependent prefactors accounting for an RG- 
scaling that leaves the full action invariant. This re- 
duces the remaining fc-dependence in r( n ' to dimension- 
less ratios x = p 2 /k 2 . Hence, the interpolation between 
the physical infrared regime (0J and the regularised be- 
haviour for p 2 k 2 is fc-independent. This implies a 
fixed point behaviour. So far, fixed point investigations 
within the exact renormalisation group have been put 
forward within derivative expansions. Here we provide a 
new truncation scheme for accessing the full momentum 
dependence of Green functions in the infrared limit. 

Taking advantage of their simplified structure in 
the physical regime (QJ, we introduce the following 



parametrisation for the gluon and ghost two point func- 
tion 

T i%(p 2 ) = Z A ■ Z A (x) ■ p 2 • n(p) • 11 + longitudinal , 
r i 2 c(P 2 ) = zc-Z c (x)-p 2 -K, (6) 
where 

Z(x) = x K (l + SZ{x)) . 

Here H^{p) = 5^ - p^/p 2 , ll ah = 8 ab and Z, 8Z, k 
without subscripts stand for Zaic-i SZa/ci k,a/c- The 
longitudinal part of the gluon two point function in © 
is not given here as it does not contribute to the Landau 
gauge flows considered below. The momentum depen- 
dence of (|fj| in the asymptotic regimes x 3> 1 and i<1 
fixes the functions SZ in both regimes: the prefactor x K 
encodes the leading behaviour for the momentum regime 
(@J where l«i<C Aq CD /fc 2 , hence SZ tends to zero. 
The fc-independence of (jBJ at k = implies dtn = and 
dtz = 2k z. In the limit of small x, we recover the trivial 
momentum behaviour T^\p 2 ) = m 2 + z\p 2 + 0(p 3 ) as 
a consequence of the cut-off. Accordingly the functions 
SZ have the limits 

SZ(x — > oo) — > , 
SZ(x^O) -» -l + px- {1+K) + 0(x~ K ). (7) 

We emphasise that the limit x — > oo implies k — > as 
we demand x <§; Aq CD /fc 2 . In the intermediate region 
x i=a 1 the propagators show a regulator-dependent inter- 
polation between the physical infrared behaviour and the 
trivial cut-off regime. Also 

d t 8Z{x) = -2xd x 5Z (8) 

reflects that the dependence of SZ on k stems from the in- 
terpolation property while its shape is independent of it. 

(n) 

The vertices Lj, in the present truncation are schemat- 
ically given by 



for n > 3. In the above parametrisation, the non- renor- 
malisation of the ghost-gluon vertex at k = leads to 

g 2 1 

k a = -2k c , a s = - T . (10) 

47T ZA Zq 

It follows from l|10(l and dtz — 2k z that the strong cou- 
pling constant achieves a fixed point in the infrared limit, 
dtOi s = 0. In particular, the z-factors for the three- 
and four-gluon vertices are proportional to <x~J 2 z 3 / 2 and 
a s z\, respectively. Moreover it follows within the trun- 
cation 10 that no mass term for the ghost is present 
corresponding to fic — in Q. In summary, the trun- 
cation © , © with the properties , ifTUj) satisfy the RG 
properties of Landau gauge QCD as well as the trun- 
cated Slavnov- Taylor identities valid in the presence of 
the regulator. 
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We proceed with determining the strong coupling fixed 
point a s and the infrared exponents k from J3J in the 
truncations put down in © . We fix the tensor structure 
of the regulators as Rf v = R A 5^S ab and R ab = R c S a \ 
leaving the scalar part R — zp 2 r(x) at our disposal. The 
integrated flow leads to an integral equation for SZ, 



5Z A/C (x) 



x 



dx' 

fA/c(x) 



(11) 



where 
f A {x) 
fc(x) 



3 Jo 7-i urc(u) + Z c {u) 



dyy dt{l-t 2 f' 2 



i (y Tc{y) 



T A {y) 



u \urA{u) + Za(u) rc{u) + Zc(u) 

and u = x + y + 2t^/xy. We have also introduced the 
abbreviation 2T = k 2 G(dtR)G, which reads explicitly 



nv) 



1 Kr(y) — yr'(y) 
y (r(y) + Z{y)) 2 



In the integration over k was traded for an integra- 
tion over x, using z = z x~ k (Aq CD /p 2 ) k with d t z = 0. 
The purely gluonic diagrams in @ are suppressed by pos- 
itive powers of p 2 /Aq CD in the fixed point regime and do 
not contribute to Ifllfl . 

For the sharp cut-off the kernel T(x) simplifies and 
becomes T(x) = S(x — l)/(zZ(l)). Then the t and y 
integrals in Ullfl can be performed analytically leading to 
a lengthy sum of hyper-geometric functions. For generic 
smooth cut-off, k and a s are computed from Ullfl in the 
limit x — > as follows. On the right hand side of (|11|) . the 
integrands /a/c are suppressed for small momenta x and 
y due to Ta/c- This implies that SZ can only contribute 
for sufficiently large arguments where it tends to zero due 
to 0. Therefore we can safely neglect SZ under the mo- 
mentum trace on the right hand side of Ijlll) in a first step. 
Then, we proceed to fully consistent solutions SZ of (jl 1|> 
by iteration. On the left hand side of the gluon equation 
the term 8Za(x — > 0) is determined by (7J. In the ghost 
equation, we have to carefully identify the terms. For 
Kc > we read-off from that SZ c {x — > 0) diverges. 
However, T { ^ c (p 2 ) = z c x KC (l + SZ c (x)) has a Taylor ex- 
pansion about x — due to the IR regularisation. Hence 
the coefficient of the constant term in SZq is —1. 

The integral equation Ullfl can be solved numerically 
for k and a s . We have restricted ourselves to the domain 
kc £ (|, 1) For general regulator r(x) we obtain 



k c = 0.59535. 



a s = 2.9717. 



(121 



which agrees with the analytic DS result |6j. This might 
seem surprising given the completely different representa- 
tions. However, the equivalence can be proven as follows: 



in the derivation of (|12fl we employed the truncation © 
and 10 with the additional constraint SZ w under the 
momentum trace in Integrating the flow © in this 
setting from to A; leads to 

z A p 2 x KA SZ A = Tr[G5 (3) G5 (3) - iGS* (4) G]J, 

z cP 2 x KC SZ c = Tr[G5 (3) G5 (3) ]J, (13) 

for the gluon and ghost two point function. Eq. 113fl al- 
ready implies an UV renormalisation at p 2 — k 2 . This is 
reflected by (|llfl being finite. Within the present trunca- 
tion the renormalised DS equations used in Q are equiv- 
alent to (fl3fl . Beyond this truncation the DS equation 
and the flow equation differ. Apart from additional t- 
derivative terms, all vertices in (|13fl get dressed as dis- 
tinguished from the DS equation, where both bare and 
dressed vertices are present. Consequently, the flow equa- 
tion © carries the correct RG properties if the trunca- 
tion for Tfc does. In particular, equation l|13fl with the 
truncation |jfj| and @ allows for a fully self-consistent 
renormalisation. Both the left and right hand sides of 
(f!3fl transform the same way under RG scalings, includ- 
ing the purely gluonic terms dropped in (|llfl . Indeed, 
it is the latter fact which makes their neglection fully 
consistent. The RG properties are most obvious for reg- 
ulators proportional to the wave function renormalisa- 
tions z, R ~ z R. Then, the propagators can be written 
as G = z~ 1 Gq where Go does not involve further wave 
function renormalisation factors. Also, the vertex func- 
tions involve the appropriate multiplicative powers of z 
following from IflUf) and 0. 

We proceed by iterating the above solution, also taking 
into account SZ ^ on the right hand side of the flow. 
Then pifl can be used for a numerical iterative determi- 
nation of SZ and k, a s . We have computed a s and k for 
various classes of regulators. The exponent kc takes val- 
ues in an interval bounded by the values (f!2fl from above 
and the iterated sharp cut-off result kc — 0.539 from 
below. However, the values l(T2^l are singled out by their 
physical relevance: k and a s depend on the regulator. 
This reflects the fact that the truncation and the 

choice of R define the approximation to the full problem. 
We expect that global extrema of K(R),a s (R) represent 
an optimised approximation based on © , J^J , see [|| . In 
general it is very difficult to translate this idea into prop- 
erties of optimal regulators or even in statements about 
the values of optimal observables. Here we present a new 
idea how to extract this information: let us assume that 
we have identified a regulator Rq which is an extremum 
in the space of regulators in the sense that 



Tr (SR 



5R> 



r 







(14) 



R=R„ 



for all SR. Eq. Ifl4l) implies that k(Ro), a s (Ro) are 
extrema, reflecting an optimisation of physical quanti- 
ties. Note that the z n (Ro) by themselves are not physi- 
cal quantities but they can be extrema up to RG scal- 
ings. However, Ifl4fl is more restrictive. For n = 2, 
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HI 41 is the additional demand, that SZ(x) is an extremal 
curve. Hence the question arises whether solutions Rq 
to 114J1 exist at all: the fixed point assumption underly- 
ing the whole analysis is valid for regulators R that only 
lead to a suppression of infrared modes. This implies a 
monotonous interpolation between the physical infrared 
regime i; 2 « p 2 « ^qcd an d the infrared regularised 
regime p 2 <C k 2 . Then the extremal curves exist and are 
given by vanishing SZ |15|. 

The above considerations allow us to identify the opti- 
mal values for k,q (R) and a s (R) without further compu- 
tation. The correction terms proportional to 5Z under 
the momentum trace are negligible near optimal regula- 
tors and we are left with the non-iterated equation 
with an optimal regulator Rq. However, we have already 
shown that the non-iterated solutions are independent of 
the regulator and we conclude that kc(Ro), cts(Ro) are 
given by (jl2|l . We are also in a position to device opti- 
mal regulators. For these regulators the effects of 8Z on 
the right hand side of (|llfl have to disappear. They need 
not to satisfy (|14fl which only applies to a small sub-set 
of optimal regulators. In Fig. 1 we show 5Z^/ C for the 
class of regulators 

r c (x) = z c 0- 1 (l-x), r A (x)= 1 * A (15) 

x(l + x) 

The 7-dependence of 5Za is negligible on the scale of 
Fig. 1. The ghost regulator in lfT5l) is a sharp cut-off 
while the gluon regulator is a rather soft mass-like cut-off. 
In the limit 7 — > 00 the solution of Qllfl approaches the 
optimal values (|12|) . One verifies that all contributions 
from SZ to the right hand side of Ijllfl tend to zero for 
7 — > 00 and we are left with the non-iterated result (|12[) . 
Hence 1150 with 7 — > 00 belongs to the class of optimal 
regulators. It is not in the extremal subset as only the 
ghost propagator satisfies (|14|) with SZc(x 7^ 0) = 0. 
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Figure 1: 5Z for regulators with 1151 

In summary, we have shown how the non-trivial mo- 
mentum dependence of Green functions in the infrared 
sector is extracted within the exact renormalisation 
group. To that end we have developed new fixed point 
and optimisation techniques applicable to this situation. 
This reasoning has been applied to the infrared behaviour 
of ghost and gluon propagators in Landau gauge QCD. 
Our results provide further evidence for the Kugo-Ojima 
scenario of confinement: the ghost propagator is infrared 
enhanced while the gluon propagator develops a mass 
gap. In a truncation with dressed propagators, we have 
recovered previous results from Dyson-Schwinger equa- 
tions in a very simple manner. This is noteworthy also 
in view of the qualitative differences between these ap- 
proaches. Our method allows for a straightforward im- 
provement of the truncation by adding vertex corrections 
and further terms with the correct RG properties. The 
same applies to dynamical quarks. In both extensions the 
correct RG scaling as well as the inherent finiteness of the 
integrated flow are most important. Further results and 
details will be presented elsewhere. 
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